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Abstract 

We determine two-particle scattering phase shifts and mixing angles for quantum theories 
defined with lattice regularization. The method is suitable for any nonrelativistic effective 
theory of point particles on the lattice. In the center-of-mass frame of the two-particle 
system we impose a hard spherical wall at some fixed large radius. For channels without 
partial-wave mixing the partial-wave phase shifts are determined from the energies of the 
nearly-spherical standing waves. For channels with partial-wave mixing further information 
is extracted by decomposing the standing wave at the wall boundary into spherical harmon- 
ics, and we solve coupled-channels equations to extract the phase shifts and mixing angles. 
The method is illustrated and tested by computing phase shifts and mixing angles on the 
lattice for spin-1/2 particles with an attractive Gaussian potential containing both central 
and tensor force parts. 
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I. INTRODUCTION 



There have been several recent studies on the subject of lattice simulations for low-energy 



nuclear physics using effective interactions 



iii,y,y,y,B,y,fly,y, 



10 



11 



12 



13|, 



14| . Similar 

lattice effective field theory techniques have been used to study cold atomic Fermi systems 
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13, 
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2l|. 



in the limit of short-range interactions and large scattering length 
In nearly all cases the connection between lattice interactions and physical observables is 
made using some variant of Liischer's result 
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23l . |2J] relating the energy levels of two- 



body states in a finite-volume cubic box with periodic boundaries to the infinite-volume 
scattering matrix. 

Liischer's method has been extended in a number of different ways. Several studies 
have looked at asymmetric boxes 
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26| . while another considered small volumes where 



the lattice length L is smaller than the scattering length 27|. There have also been studies 



of moving frames 



28|, 



29| . winding of the interaction around the periodic boundary 30|, 
modifications at finite lattice spacing , and techniques to distinguish shallow bound states 
from scattering states using Levinson's theorem [3l]]. A very recent study derived finite 
volume formulas for systems of n bosons with short-range repulsive interactions 32|. 

While Liischer's method has been very useful at low momenta, there is currently no 
technique which is able to determine phase shifts on the lattice at higher energies and higher 
orbital angular momenta. There is also no technique which can accurately measure spin- 
orbit coupling and partial-wave mixing on the lattice. The physics of spin-orbit coupling and 
partial-wave mixing is difficult to extract using Liischer's method due to artifacts generated 
by the periodic cubic boundary. An example of this problem is shown in the summary 
and discussion section of this paper. Any method overcoming these theoretical problems 
would probably not be numerically practical when applied to lattice simulations of quantum 
chromodynamics. This is because each individual hadron must be constructed as a bound 
state of quark and gluon fields, and the Monte Carlo signal for hadron-hadron scattering 
states well above threshold would be very weak. However for an effective lattice theory 
of fundamental point particles, such a technique should be numerically viable since the full 
two-particle spectrum is relatively easy to compute. 

In this paper we discuss a simple method which directly measures two-particle phase shifts 
and mixing angles on the lattice. In the center-of-mass frame of the two-particle system 
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we impose a hard spherical wall boundary condition at some large radius. Phase shifts 
and mixing angles are determined from properties of the nearly-spherical standing waves 
produced by the wall boundary. The method is suitable for any nonrelativistic effective 
theory of point particles on the lattice with or without spin. We test the method on the 
lattice using an attractive Gaussian potential for spin-1/2 particles containing both central 
and tensor force parts. 

The organization of the paper is as follows. We first discuss representations of the cubic 
rotational group. We then introduce the test potential and solve the S = channels and 
uncoupled S = 1 channels up to J = 4 using the spherical wall method. We then discuss 
coupled partial waves with a common nodal constraint and the coupled equations needed to 
solve for the L = J — 1 and J + 1 phase shifts and mixing angle for total angular momentum 
J. Finally we solve the coupled S = 1 channels on the lattice up to J = 4 and discuss 
applications of the technique. 



II. CUBIC ROTATIONAL GROUP 



Our choice of lattice regularization reduces the SO (3) rotational symmetry of continuous 
space to the cubic rotational group SO (3, Z). This consists of 24 group elements generated 
by products of 7r/2 rotations about the x, y, z axes. Since SO (3, Z) is discrete, we cannot 
define angular momentum operators J^, Jy, Jz in the usual sense. However if R-^ {t^/'^) is 
the group element for a 7r/2 rotation about the z axis, we can use the SO (3) relation 

i?,(7r/2) = exp [-z^J,] (1) 

to define Jz- The eigenvalues of Jz are integers modulo 4. and Jy can be defined in the 
same manner. 

There are five irreducible representations of the cubic rotational group. These are usually 
written as Ai, Ti, T2, and A2. Some of their properties and examples in terms of spherical 
harmonics YL^L^{9,(f)) are listed in Table [H 

The 2 J + 1 elements of the total angular momentum J representation of SO (3) break up 
into smaller pieces consisting of the five irreducible representations. Examples for J < 7 
are shown in Table HTl 33 1. 

In this analysis we consider the scattering of two identical particles and assume that 
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TABLE I: Irreducible S0(3,Z) representations 



Representation 


Jz 


Example 


Ai 


0mod4 




Ti 


0,1, 3 mod 4 


{yi,o,ii,i,ii,-i} 


E 


0, 2 mod 4 


\v ^^2,-2+^2,2 1 
1^2,0, ^ ] 


T2 


1,2, 3 mod 4 




A-2 


2 mod 1 





TABLE II: S0(3, Z) decompositions for J < 7 



S0(3) 


S0(3,Z) 


J = 


^1 


J = 1 


Ti 


J = 2 


E®T2 


J = 3 


Ti e Ta e ^2 


J = 4 


Ai e Ti ® E © T2 


J = 5 


Ti®Ti®E®T2 


J = 6 


Ai®Ti®E®T2®T2®A2 


J = 7 


Ti®Ti®E®T2®T2®A2 



the interactions are parity conserving. Therefore in continuous space the two particles 
share the same group representation for intrinsic spin and the same group representation 
for all other internal quantum symmetries. As a result each two-particle state must be 
symmetric or antisymmetric under parity and symmetric or antisymmetric with respect 
to internal symmetry quantum numbers. Once these are fixed the overall Fermi or Bose 
statistics determines whether the total intrinsic spin combination must be symmetric or 
antisymmetric. 

There are well-known problems associated with massless fcrmions on the lattice. How- 
ever for nonrelativistic particles one can easily find a lattice Hamiltonian or Euclidean action 
which maintains exact parity invariance and all internal symmetries. Therefore the sym- 
metry or antisymmetry of intrinsic spin on the lattice is unambiguous and the same as in 
continuous space. For two spin-1/2 particles where the total intrinsic spin is S — Q or 
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S — 1, this is enough to specify the intrinsic spin representation completely. If the intrinsic 
spin is symmetric then the S0(3, Z) representation is Ai. If the intrinsic spin is antisym- 
metric then the SO (3, Z) representation is T^. These are in one-to-one correspondence with 
the representations S* = and S* = 1 in continuous space, respectively, and so there is no 
confusion in borrowing the continuous space names for the two cases. For particles with 
higher intrinsic spin, there is in general some unphysical mixing on the lattice among even 
values of S and among odd values of S. The same mixing on the lattice occurs among even 
values and odd values of orbital angular momentum L regardless of the intrinsic spin. 



III. TEST POTENTIAL 



There is an endless variety of different interactions one can study on the lattice. For our 
analysis we choose a simple pedagogical example for which the continuum limit is relatively 
easy to compute. We choose a bounded short-range potential producing both a central 
force and tensor force. This yields all of the essential comphcations of spin-orbit couphng 
with partial-wave mixing while leaving out additional issues of ultraviolet divergences and 
singular interactions. We take a system of identical spin- 1/2 particles with mass m and a 
spin-dependent potential, 

V{f) = {l + ^ [3 (f ■ aO (r ■ a,) - a, ■ a^lj exp (-^^) • (2) 

The range of the potential is set by the parameter Rq. We leave the internal symmetry 
group unspecified and consider all possible combinations of symmetry or antisymmetry for 
parity and intrinsic spin. 
The tensor operator 

'S'i2(r) = 3 (f • ai) (f • CT2) - (Ti • ^2 (3) 

is of considerable general interest in physics. It is produced both by one-pion exchange in 
low-energy nuclear physics as well as by magnetic and electric dipole interactions in atoms 
and molecules. The parameters C, Rq, and m we choose are motivated by low-energy 
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nuclear physics, 



C = -2 MeV, (4) 
i?o = 2 X 10-2 MeV-\ (5) 
m = 938.92 MeV. (6) 

This produces a very shallow bound state in the ^S{D)i channel with energy —0.155 MeV. 
The parenthesis in our spectroscopic notation indicates a mixture of two partial waves. For 
example ^S{D)i indicates a mixture of S and D partial waves. The ordering of the partial 
waves, S{D) rather than D{S), indicates that this state becomes a pure 5* wave when the 
tensor part of the interaction is continuously dialed down to zero. This identification is 
unambiguous in finite volume systems where the energy levels are discrete. 

For total intrinsic spin S = the tensor operator Su vanishes. And so for this case it 
suffices to consider just the central part of the potential, 

1 r2 

i 

for particles with no intrinsic spin at all. We start first with this simplest case. Afterwards 
we consider S = 1 for uncoupled channels and then finally 5 = 1 in coupled channels. 

IV. INTRINSIC SPIN 5 = 

We measure phase shifts by imposing a hard spherical wall boundary on the relative 
separation between the two particles at some chosen radius -Rwaii- Viewed in the center-of- 
mass frame we solve the Schrodinger equation for spherical standing waves which vanish at 
r = -Rwaii- This is sketched in Fig. [1] We are of course assuming that the original lattice 
system is large enough to hold a sphere of radius -Rwaii- The reason for this spherical wall 
is to remove copies of the two-particle interactions due to the periodic boundaries of the 



Vo{n = Cexp{--^], (7) 



lattice. It is this feature that distinguishes our method from Liischer's method \24i\ . 

For values of r beyond the range of the interaction, the spherical standing wave in con- 
tinuous space can be decomposed as a superposition of products of spherical harmonics and 
spherical Bessel functions, 

[cos Sl X jiikr) - sinSi x yL{kr)] Yl,lA^A), (8) 
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FIG. 1: Spherical wall imposed in the center-of-mass frame. 



where the center-of-mass energy of the spherical wave is 



E = 2x 



2m 



m 



(9) 



and the phase shift for partial wave L is 5^. Therefore we know k from the energy E, and 
the phase shift 6l is determined by setting the wavefunction ^ equal to zero at the wall 
boundary, 

cos Sl X iL{kR^^\\) = sinSi x yiikR^^n), (10) 

Sl = tan — — . 11) 

On the lattice there is some ambiguity on the precise value of -Rwaii since the components 
of r must be integer multiples of the lattice spacing. We resolve this ambiguity by fine- 
tuning the value of -Rwaii for each standing wave so that 6l equals zero when the particles 
are non-interacting. This is illustrated in the following discussion of the lattice calculation. 



V. LATTICE CALCULATION FOR S = 



Since our method is intended as a tool for simulations of few- and many-body systems 
using lattice effective field theory, we present the analysis using the standard formalism of 
lattice field theory developed in the literature. For the free part of the lattice action we 
use the same lattice action defined in |l^ with spatial lattice spacing a = (100 MeV)~^ 
and temporal lattice spacing at = (70 MeV)"^. We define at as the ratio between lattice 
spacings, at = at/ a. Throughout we use dimensionless parameters and operators, which 
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correspond with physical values multiplied by the appropriate power of a. Final results are 
presented in physical units with the corresponding unit stated explicitly. 

Since the temporal lattice spacing is nonzero we work with transfer matrices rather than 
the Hamiltonian directly. In simple terms the transfer matrix is just the exponential of the 
Hamiltonian exp(— iJAt), where At equals one temporal lattice spacing. More precisely 
the free-particle transfer matrix is defined as 



free 



;: exp (-iffreca*) : , 



(12) 



where the :: symbols indicate normal ordering. We use the O(a^) -improved free lattice 
Hamiltonian, 

49 



free 



12m 



n j=^,^ 



+ 



3 

i^HHH [a]{n)a,{n + 2l)+a]{n)a,{n~2l) 
5Z [a]{n)aj{n + 31) + a] (n) a,- (n - 3/) 



H jr-=t,J, (=1,2,3 



(13) 

n j=T,i (=1,2,3 '" 

The vector n denotes integer- valued coordinate vectors on a spatial three-dimensional lattice, 
and / = i, 2, 3 are lattice unit vectors in each of the spatial directions. 

Since the potential energy is finite we can take the lattice potential at n to agree with 
the continuum potential at r = na. For = we can omit the tensor part of the potential, 
and so the transfer matrix is 



M = : exp 



-//free"* - y 5Z - n2)p''''%n,)p''''''{ 



n2 



ni,n2 



(14) 



where p"-^'"-{n) is the particle density operator 



j=T,i 



n)aj[n). 



(15) 



3^. In the 



We calculate eigenvalues of the transfer matrix using the Lanczos method 
transfer matrix formalism eigenvalues of the transfer matrix are interpreted as exponentials 
of the energy, 

M|^) = A|^) = e"^"* 1^). (16) 
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FIG. 2: Free particle spectrum for standing waves with radius -Rwaii = 10 + e lattice units. 

In Fig. [2] we show the free two-particle energy spectrum in the center-of-mass frame. We 
have set the spherical wall so that the amplitude for two-particle separation greater than 
10 lattice units is strongly suppressed by a large potential energy due to the spherical wall. 
We can write this as 



Vo{ni - ns) Vo{ni - ng) + Kvaii x 



|ni -nsl - (10 + e)), 



(17) 



where 6 is the unit step function and e is a small positive number. This choice for the 
spherical wall radius is large enough to probe wavelengths much larger than the lattice 
spacing. The breaking of rotational invariance due to lattice regularization can be seen in 
the small splitting of the different SO (3, Z) representations comprising each orbital angular 
momentum multiplet. The elements of each S0(3, Z) representation however remain exactly 
degenerate. 

In Fig. [3] we show the interacting two-particle energy spectrum for S = for the same 
wall radius. From the lattice data we find that the interaction is attractive in each of the 
S = channels. This is as one might expect since Vq is negative definite. To illustrate the 
calculation of phase shifts, we take for example the 1^5*0 state. For the free particle system 
the 1^5*0 state has energy 0.9280 MeV, which corresponds with a continuum momentum of 



fcfree = 29.52 MeV. 



(18) 



> 



5 = 0, even parity 
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L, mod 4 
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5 = 0, odd parity 
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FIG. 3: Interacting spectrum for 5 = standing waves with i?waii = 10 + e lattice units. 



The 1 5*0 state has exactly one extremum and so the free wavefunction vanishes when 

jo(/Cfrcc^wall) = 0, (19) 

i?waii = = 0-1064 MeV-^ (20) 

This corresponds with 10.64 lattice units, and is consistent with our wall starting at lattice 
distances greater than 10 lattice units. For the interacting system the 1^5*0 state has energy 
0.6445 MeV. This corresponds with a continuum momentum of 

k = 24.60 MeV. (21) 

Using (fTT!) we find a ^5*0 phase shift at A; = 24.60 MeV equal to 

^("'^S'o) = tan^ 

We proceed in this manner for all of the S = phase shifts. In cases where the spin 
multiplet is slightly non-degenerate due to the lattice we use the average energy over all 
elements of the multiplet. 

In Fig. Hlwe compare lattice results and exact continuum results for the ^5*0, ^Pi, ^D2, 
^F^, and phase shifts as a function of the center-of-mass momentum pcu- For the lattice 
results we use -Rwaii = 10 + e, 9 + e, and 8 + e lattice units. In order of increasing momentum, 

10 



UkR 



wall J 



VoikR 



wall J 



30.0°. 



(22) 
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FIG. 4: 5 = phase shifts for J < 4. In order of increasing momentum the lattice data corresponds 
with the first radial excitation for i?waii = 10 + e,9 + e, and 8 + e; second radial excitation for 
-Rwaii = 10 + e, 9 + e, and 8 + e; and so on. 



the lattice data corresponds with the first radial excitation for -Rwaii = 10 + e, 9 + e, and 
8 + e; second radial excitation of -Rwaii = 10 + e, 9 + e, and 8 + e; and so on. The continuum 
results are calculated by solving the Lippmann-Schwinger equation in momentum space. 
The lattice results are within a few percent of exact results for momenta below 80 MeV. 
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The error increases to 10% or 15% for momenta near 120 MeV. This is as good as can be 
expected without fine-tuning of the lattice action. The lattice spacing a = (100 MeV)~^ 
corresponds with a momentum cutoff equal to vr/a = 314 MeV. The quality of the results 
at higher spin is remarkable considering that can only be determined modulo 4 on the 
lattice using ([T]). 

VI. LATTICE CALCULATION FOR 5 = 1 

For intrinsic spin S = 1 the tensor interaction in V{r) must be included. We define the 
spin densities, 

pf'^in) = «:' (^) ^ = 1' 2> 3. (23) 

The transfer matrix has the form 



M = : exp 



ni,n2 



at 



2YI Yl M^i-n2)Tu>{ni-n2)pf'''ini)pf-[n2) 




a' .a / 



2R. 

ni,n2 l,l'=l,2,3 



(24) 



where 



Twin) = 3nini> - \nf 5w. (25) 

In Fig. [5] we show the interacting two-particle energy spectrum for 5* = 1 for wall radius 
-Rwaii = 10 + e lattice units. 

The total angular momentum J multiplets are deduced from the approximate degeneracy 
of S0(3, Z) representations comprising the multiplet decompositions in Table HIl In some 
cases the accidental degeneracy of different J multiplets makes this process difficult. For 
example the 1^F{P)2 and l^F{H)i multiplets are nearly degenerate for -Rwaii = 10 + e, as 
can be seen in Fig. [51 In such cases further information can be extracted by calculating 
the inner product of the standing wave with spherical harmonics Yl l_,{9, (p). In the infinite 
volume limit the l^S{D)i bound state has an energy of —0.170 MeV, in good agreement 
with the exact result —0.155 MeV. 

For S* = 1 we start with the uncoupled channels. Once the "^^^^Lj multiplets are 
identified the process is exactly the same as for the S* = case. In Fig. [6] we compare 
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FIG. 5: Interacting spectrum for 5 = 1 standing waves with Rwaii = 10 + e lattice units. 



lattice results and exact continuum results for the ^Pq, ^Pi, ^D2, ^F^, and phase shifts. 
For the lattice results we again use -Rwaii = 10 + e, 9 + e, and 8 + e lattice units. 

The tensor interaction produces a strong repulsion in the ^Pq channel, enough to overcome 
the attraction from the central force interaction. As in the S = case, the lattice results 
for the uncoupled S = 1 channels are within a few percent of exact results for momenta 
below 80 MeV and within 10% or 15% for momenta near 120 MeV. 



VII. PARTIAL WAVE MIXING FOR S = 1 AT ASYMPTOTICALLY LARGE RA- 
DIUS 



In the next few sections we consider partial-wave mixirig between the partial waves L = 
J — 1 and L = J + 1 using the Stapp parameterization [35] . It is convenient to use the 
two- component notation, 

^ Rj-i{r) 



(26) 
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FIG. 6: Uncoupled S = 1 phase shifts for J < 4. The lattice data corresponds with 
10 + e,9 + e, and 8 + e. 



for the radial part of the wavefunction in continuous space. The full expression is 

j-i 1 

Lz=-{J-1) S,=-l 
J+1 1 

+ i?j+i(r) Yl E yj+i,LAOA){J+l,W,l,S,\J,J,)x\l,S,), 



14 



where {L, L^; S, Sz \J,Jz) denotes the usual Clebsch- Gordon coefficient for adding orbital 
angular momenta and intrinsic spin. Using this shorthand notation the S-matrix can be 
parametrized as a 2 x 2 matrix of the form 







COS 2ej i sin 2ej 
ism2ej cos 25 j 



e'^''+^ 

Let the incoming wave at asymptotically large radius r be 

^ in 





Q 




C 


2ikr 


e*'^/^ 




D 



Then the outgoing wave at large r is [35|, [Sg, [37 1 



out 



J{kr-jTT/2) 

2ikr 
2ikr 



Jtt/2 g 







-i7r/2 



c 

D 



Q^in/2^2iSj_i ^.Qg 2ej + iL)e^''/2g*(5,7-i+5j+i) gjj^ 2ej 



It is convenient to define another unitary operator 



W 







cose J —tsmej 
—i sine J cosej 



so that 



cose J zsmej 
i sine J cosej 







AS ] 











D«0J + 1 



cose J ismej 
i sine J cose J 



These satisfy the identities 



WW^ = W^W = 1, 
W*W^ = S. 

To construct a real-valued standing wave, we let the asymptotic incoming wave be 

2ikr 



(28) 



(29) 



(30) 



(31) 



(32) 
(33) 

(34) 
(35) 



g-i7r/2 Q 




c 


e'^''^ 


w 






D 



(36) 
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for real numbers C and D. Then the outgoing wave is 

gi(A;r-j7r/2) 



2ikr 



Q 





^ = 2Re 
1 

kr 



^i{kr-J'K/2) 

2ikr 





C 


i(kr-JiT/2) 


sw 








L> 




2ikr 


is then 










gi7r/2 







C 






iTr/2 


W* 




1 


e- 




D 





Q 




C 


e-^'^/^ 











(37) 



C sin (/cr — "^tt + cosej + D cos {k:r — ^^7i + sinej 
D sin (/cr — -^vr + cos ej + C cos (/cr — -^tt + 5j+i) sin £ j 



(38) 



Wc now impose a hard spherical wall boundary at r = i?waii- Both partial waves vanish 
at the wall, and we define the angle A so that 



Then 



and 



C 
D 



- A = A;i?wall Ti—TT. 



tan (—A + 5j-i) — — tan£j 



^ tan (-A - TT + 5j+i) = ^ tan (- A + Sj+i) = - tan ej. 



(39) 

(40) 
(41) 



In general there are two solutions for A per angular interval vr. In the asymptotic case where 
-Rwaii ^ the spacing between energy levels becomes infinitesimal. We can therefore 
choose two independent solutions A^ and A^^ with nearly the same values of k, and we 
neglect the difference in k in the following steps. 
The solutions A^ and A^^ satisfy 



tan (- A^'-^^ + tan (- A^'-^^ + 5j+i) = tan^ Sj. 
The standing wave now looks hke 



kr 

C_ 
kr 



(42) 



C sin (/cr - + 5j_i) cos sj - C^^p^^TTT ~ + ^J-i) cos e j 

-C sin {kr - ^TT + 8j+^) l^lZ^+t'+i) ^-^ + ^ (^^ - + sin £ j 
, , sin \kr - ^71 + A] 

cos(— A+flj_i) L 2 J 



sin(- 



^^sin [kr-^n + A]_ 



(43) 
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There are two linearly independent real-valued standing wave solutions. We write them as 



kr 



AY_[ sin [kr - ^vr + A^-^^) 
AY^[ sin {kr - ^Tx + A^'^^) 



Then 



AYJi tan e j 



-A 



,^,jsin(-A^'^^ + 5j+i) 



^+^os(-A^'^^ + 5j_i: 



(44) 



(45) 



Only three of the constraints in fH2|) and (145!) are needed to determine the parameters 
and ej. Therefore one constraint, for example 

sin (-A^^ + 5j+i) 



J^j_x tanej 



(46) 



cos (-A^^ + (5j_i) ' 

must be redundant. This makes sense in light of the symmetry associated with the inter 
change of the two solutions A"^ and A"^^. The orthogonality of and implies 



A^ 



A^^ 
A^^ ' 



(47) 
(48) 



From this one can derive (HEl) from fH2l) and (HSl). 



VIII. EXPANSION FOR SMALL ej 

For our test potential we have chosen the tensor force to be very strong, and so the mixing 
angles are quite large. In several real physical systems however, for example nucleon-nucleon 
scattering at low momenta, the mixing angles are quite small. For small Ej it is useful to 
expand around ej = 0. We choose A-'^ to be the solution that equals at zero mixing, 
and A^^ to be the solution that equals at zero mixing. Since 

tan (-A^ + 5j„i) tan (-A^ + 5j+i) = tan^ £j, (49) 

the first correction to A^ and A^-'^ comes at order 0(ej), 

A^ = 5^_i + cV, + 0(e^^), (50) 
A^^ = 5j+i + c^V, + 0(e^^). (51) 

We find 

- c'^j tan (5j+i - 5j_i) = + o(£^^), (52) 
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Similarly 



Since 



we have 



Up to 0{e]) then 



or 



Similarly we have 



and so 



1 



ji 



tan(5j+i - 



J J sin (-A^ + 6j+i) 



sin - 6j-i 



0(4) 



sin (A^^- AO = -4^^ + 0(4), 



J+i 



j+i 



sin (A^^- A^ 



^j^itan^j — — 



sm 



4^^ — —A^^ 



II 



cos (-A^^ + 



J 



£jcos(-(5j+i + 



sin(5j+i - 6j-i) 



or 



Up to order 0(£:j) we find 



= 4^ sin (A^^- AO +0(4). 



This is consistent with fl59l) given the orthogonality condition, 
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IX. HARD SPHERICAL WALL AT NON-ASYMPTOTIC RADIUS 



The constraints in (H2|) and (H5|) hold if the spherical wall is at asymptotically large 
radius, fc-Rwaii ^ 1- In this limit the nodes of the spherical Bessel functions jj-i{kr) and 
jj+i{kr) coincide, and we determine one shift angle for both partial waves, 

J- 1 



A = kR 



wall 



-vr. 



(65) 



However for numerical calculations it is more convenient to work with smaller values of 
kRwaii- In this case the coincidence of the nodes for L = J — 1 and L = J + 1 at the wall 
boundary does not automatically imply the coincidence of nodes for kr ^ 1. Nevertheless 
we can still work with the asymptotic form for kr ^ 1 if the node remains at i?waii but the 
spherical wall is removed. 



1 

kr 



C sin {kr — "^vr + cos Sj + D cos {kr — "^tt + sin Sj 

D sin {kr — "^tt + 5 j+i) cos 6j + C cos {kr — ^^ix + sin sj 



(66) 



Instead of one angle A, we define two angles Aj_i and Aj+i so that the J — 1 partial wave 
vanishes when 

- Aj_i = kr- ^^^TT, (67) 



and the J + 1 partial wave vanishes when 

- Aj+i = kr 

Then 

C 



J+1 



-TT. 



and 



— tan(-Aj_i + 5j^i] 



— tan(-Aj+i + bj+i, 



— tan Ej 



tanej. 



(68) 



(69) 



(70) 



In general for a given /c, there are two pairs of solutions (Aj_i, Aj_|_i) per angular interval 
TT corresponding with different positions of the nodes. However instead of taking k as the 
independent variable, we consider the location of the spherical wall radius as the independent 
variable. We denote the two solutions as (Aj_j^, Aj,^^, k^^ and (Aj^^^, Aj{^]^, A;^^) . We have 



tan (-Aj_i + 6j-i) tan (-Aj^^ + 5j+i) = tan^ej. 



(71) 
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at momentum , and 



tan [—Aj^_^ + tan (— Aj^^^ + = tan^ej, 



(72) 



at momentum k^^ . 

The standing wave now looks like 



1 

kr 

_C 

kr 



C sin {kr — "^tt + cos ej + D cos (A;r — "^tt + sin €j 
D sin (A;r — ^^vr + cos ej + C cos [kr — ^^tt + sin ej 



cos £j 



sm 



[kr - ^vr + Aj_i] 



sin(-Aj+i+(5j+i) 

We write these as 



cos(-Aj_i+<5j_i) 

sin^, _ _ J±l^ ^ ^^^^j 



OC 



sin [k^r - ^tt + A5_i) 



A j+i sin 



J+i 



TT 



and 



Ay_^ sin (fc^V - + A5^_i) 



y45{^i sin {k^^r — 



J+i 



(73) 



(74) 



(75) 



So the generalization of equations ( l42l) and ( l45l) to the non-asymptotic case is 
tan + (5j-i) tan {-^j+i + 5j+i) = tan^ej, 

tan (-Ay_i + 5j_i) tan (-A^+i + 5j+i) = tan^ej, 



Aj_^ tanej = —Aj 



J+i 



Aj_-^ tan £^ J" — Ajj^Y 



II 



sin (-A^^^ + (5j+i) 
cos (-A5_i + 5j_i) 

sin (-Ay+i + (5j+i) 



cos 



(-Ay_i + 5,_i)' 



(76) 
(77) 
(78) 

(79) 



We note that the phase shifts and mixing angle in fl76p and (1781) are at momentum k^ while 
the phase shifts and mixing angle in (!77|) and (1791) are at momentum k^^ . In the calculation 
one must in general interpolate between values for k = k' and k = k^^ . However the need 
for interpolation is largely eliminated by considering only close pairs of values k^ ~ k^^ in 
solving ( !76l) -( l79l) . This can be done for example by considering the (n+l)^*-radial excitation 
of L = J — 1 together with the n^'^-radial excitation of L = J + 1. In this scheme we use 



tan (-A:^_i + 6j-i{k^)) tan {-A^j^^ + 6j+i{k^)) = tan^ [ej{k^)] , 
tan (-Ay_i + 6j-i{k^)) tan (-A^+i + 6j+i{k^)) ^ tan^ [£j{k^) 
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(80) 
(81) 



cos (-A^j_i + ()j_i(A;^)j 



(82) 



for the phase shifts and mixing angle ai k = k\ and 

tan (-A5_i + 5j^,{k'')) tan (-A^+i + (5j+i(fc^0) - tan^ , (83) 

tan (-Ay_i + tan (-A^+i + 5j+i(A;^0) = tan^ [ej{k'')\ , (84) 

A,_,tan [ej{k )J - + ' (85) 

for the phase shifts and mixing angle at A; = k^^ . This is the technique we use for the lattice 
results presented here. If more accuracy is required then this pair technique for k^ ^ k^^ 
can be used as a starting point to determine numerical derivatives for the phase shifts and 
mixing angle. Then the equations (175]) - fl7^ can be solved again while keeping all terms at 
order 0{k^ — k^^). 

In the limit of small mixing the expressions in ( l80i) -( l85ll reduce to 

6j_,{k^) = A^j_, + f'^^'^ + 0{e% (86) 

tan + 6j+i{k^)) 



at k = k^ and 



^Ak') = sin (Ay^, - A^,+,) + 0(6% (87) 



All 

^Ak'') = ^ Sin (Ay_, - AS_0 + 0(4). (89) 



citi k — Iv 



X. RESULTS FOR 5 = 1 COUPLED CHANNELS 

In Fig. [7] we show lattice and continuum results for ^Si, ^Di partial waves and J = 1 
mixmg ang le El. The ^P2, ^^^2 partial waves and mixing angle 82 are shown in Fig. [SI The 
■^Da, ^G^ partial waves and mixing angle £3 are shown in Fig. O The "^^4, ^Hi^ partial waves 
and mixing angle £4 are shown in Fig. (TDl As before we use -Rwaii = 10 + e, 9 + e, and 8 + e 
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20 40 60 80 100120140 20 40 60 80 100120140 

PcM (MeV) (MeV) 

El 



lattice □ 
continuum 




20 40 60 80 100 120 140 
PCM (MeV) 



FIG. 7: Coupled partial waves and mixing angle for J = 1. The pairs of points connected by 
dotted lines indicate pairs of solutions at A; = A;^ and k = k^^ . 



lattice units. The pairs of points connected by dotted hnes indicate pairs of solutions at 
k = and k = k^^ . The partial- wave ratios 

(90) 

Al ' All ' ^^^^ 

are determined by computing the inner product of the standing wave near the spherical wall 
with spherical harmonics. 

The lattice results in the coupled channels are somewhat less accurate than the S = 
results and uncoupled S = 1 results. However they are still within 20% of the exact results 
up to momentum 120 MeV. The leading error appears to come from mixing with extraneous 
channels due to broken rotational invariance on the lattice. For example the ^Si-^Di 
standing waves have a small but non- negligible amount of a ^D^ component. Nevertheless 
overall the lattice results are remarkably good for partial-wave mixing at higher total angular 
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FIG. 8: Coupled partial waves and mixing angle for J = 2. 



momentum. 



XI. SUMMARY AND DISCUSSION 



We have discussed a general technique for measuring phase shifts and mixing angles for 
two-particle scattering on the lattice. In the center-of-mass frame we impose a hard spherical 
wall at large fixed radius -Rwaii- For channels without mixing we identify total angular 
momentum J multiplets by the approximate degeneracy of SO (3, Z) representations and 
calculate phase shifts from the energies of the spherical standing waves. For channels with 
partial-wave mixing, further information is extracted by decomposing the standing wave at 
the wall boundary into spherical harmonics. The coupled-channels equations are then solved 
to extract the phase shifts and mixing angles. The method was tested by computing phase 
shifts and mixing angles for J less than or equal to 4 for an attractive Gaussian potential 
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FIG. 9: Coupled partial waves and mixing angle for J = 3. 



with both central and tensor force parts. The Gaussian envelope had a characteristic size 
of i?o = 2 X 10~2 MeV-\ and the strength was tuned to produce a shallow ^S{D)i bound 
state. At spatial lattice spacing a = (100 MeV)~^ and temporal lattice spacing at = (70 
MeV)~^ we found agreement with the exact results at the 10% — 20% level for momenta up 
to 120 MeV. 

The hard spherical wall removes extra copies of the two-particle interactions due to the 
periodic boundaries of the lattice. This appears to be important for measuring phase shifts 
and mixing angles at higher energies and higher angular momentum. For comparison we 
show in Fig. [11] the same test potential spectrum for intrinsic spin 5* = 1 on a 12 x 12 x 12 
lattice with periodic boundaries. If we compare with Fig. [5] we see that the SO (3, Z) 
representations in the various total angular momentum J multiplets have been split apart 
much further. It is easy to identify the lowest-lying standing waves but difficult to clearly 
discern much beyond this. 
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FIG. 10: Coupled partial waves and mixing angle for J = 4. 
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FIG. 11: Interacting spectrum for S = 1 standing waves for a 12 x 12 x 12 periodic lattice. 



25 



The method we have presented can be apphed directly to any nonrelativistic effective 
theory of point particles on the lattice. There is a long list of interesting few- and many- 
body systems one could consider on the lattice. For example one interesting system is that 
of atoms and molecules with long-range dipole interactions, 

1 



dipole-dipole 



(r ) OC 



47rr3 



[3 (r ■ ill) (r • /xa) - /2i ■ /I2] 



(91) 



There is interest both in magnetic dipole interactions of cold atoms 
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45 



38 



39 



40 



41 



4^ and 



46| . These spin-changing 



electric dipole interactions in cold polar molecules [43|, 
interactions are important for determining topological structures of the many-body ground 
state as well as practical issues impacting evaporative cooling in magnetic traps. 

We will discuss in some detail the application to low-energy nuclear physics in a forth- 
coming paper. One important low-energy interaction between nucleons is the exchange of 
a virtual pion, producing spin- dependent forces of the form 
2 

(ti • T2) 



^i.(r) = ( ^ 



mte 



127rr 



Si2{r) 1 + 



+ 



+ (Ti ■ (72 



(92) 



Here r are Pauli matrices in isospin space, m^r is the pion mass, /vr is the pion decay 
constant, qa is the nucleon axial charge, and Si2{r) is the tensor operator for spin-1/2 



particles defined in_ 
can be found in 




3]). Recent reviews of literature relating to chiral effective field theory 
48l |. This tensor interaction has some general similarities to that of 
the test potential considered here. However the divergent short- distance behavior of the 
one-pion exchange potential and other short-distance nucleon-nucleon interactions means 
that the lattice spacing plays an essential role in regulating ultraviolet divergences. The 
hard spherical wall method we have presented should therefore be useful in determining the 
underlying physics of a given lattice action. The singular short-distance interactions are 
likely to cause larger breaking of rotational invariance than that found for the bounded test 
potential. In such cases it may be useful to include higher-order derivative interactions 
specifically designed to cancel lattice artifacts. 
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